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Abstract 

We analyze the crystallization and collapse transition of a simple model for flexible polymer 
chains on simple cubic and face-centered cubic lattices by means of sophisticated chain-growth 
methods. In contrast to bond-fluctuation polymer models in certain parameter ranges, where these 
two conformational transitions were found to merge in the thermodynamic limit, we conclude from 
our results that the two transitions remain well-separated in the limit of infinite chain lengths. 
The reason for this qualitatively distinct behavior is presumably due to the ultrashort attractive 
interaction range in the lattice models considered here. 

PACS numbers: 05.10.-a, 87.15.Aa, 87.15.Cc 



* E-mail : Thomas . Vogel@itp . uni-leipzig . de 
'E-mail: Michacl.Bachmann@itp.uni-leipzig.de 
+E-mail: Wolfhard.Janke@itp.uni-leipzig.de 



Homepage: http : //www . physik . uni-leipzig . de/CQT . html 



1 



I. INTRODUCTION 



The analysis of conformational transitions a single polymer in solvent can experience 
is surprisingly difficult. In good solvent (or high temperatures), solvent molecules occupy 
binding sites of the polymer and, therefore, the probability of noncovalent bonds between 
attractive segments of the polymer is small. The dominating structures in this phase are 
dissolved or random coils. Approaching the critical point at the temperature, the polymer 
collapses and in a cooperative arrangement of the monomers, globular conformations are 
favorably formed. At the point, which has already been studied over many decades, the 
infinitely long polymer behaves like a Gaussian chain, i.e., the effective repulsion due to 
the volume exclusion constraint is exactly balanced by the attractive monomer- monomer 
interaction. Below the temperature, the polymer enters the globular phase, where the 
influence of the solvent is small. Globules are very compact conformations, but there is 
little internal structure, i.e., the globular phase is still entropy-dominated. For this reason, 
a further transition towards low-degenerate energetic states is expected to happen: the 
freezing or crystallization of the polymer. Since this transition can be considered as a 
liquid-solid phase separation process, it is expected to be of first order, in contrast to the 
transition, which exhibits characteristics of a second-order phase transition l|, |2J. 

The complexity of this problem appears in the quantitative description of these processes. 
From the analysis of the corresponding field theory Q| it is known that for the transition 
the upper critical dimension is d c = 3, i.e., multiplicative and additive logarithmic corrections 

n n n a 

to the Gaussian scaling are expected and, indeed, predicted by field theory [4J, [5|, la, 17[. 
However, until now neither experiments nor computer simulations could convincingly provide 
evidence for these logarithmic corrections. This not only regards analyses of different single- 



polymer models 



10, 



11 



12 



unmixing in polymer solutions 



13. 



16 



14 [jjput a so the related problem of critical mixing and 



Ini,ii8i,ii9|. 

In a remarkable recent study of a bond-fluctuation polymer model, it was shown that, 
depending on the intramolecular interaction range, collapse and freezing transition can fall 



together in the thermodynamic limit 12j, [l3| . This surprising phenomenon is, however, not 



general. For an off-lattice bead-spring polymer with FENE (finitely extensible nonlinear 
elastic) bond potential and intra-monomer Lennard- Jones interaction, for example, it could 

n 

be shown that both transitions remain well separated in the limit of infinitely long chains |14| . 
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In our study, we investigate collapse and freezing of a single homopolymer restricted to 
simple cubic (sc) and face-centered cubic (fee) lattices. We primarily focus on the freezing 
transition, where comparatively little is known as most of the analytical and computational 



studies in the past were devoted to t 



e.g., Refs. 
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16 



3, q, y, 0, 



ie controversially discussed collapse transition; see, 
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26l |. A precise statistical analysis of the 



conformational space relevant in this low-temperature transition regime is difficult as it is 
widely dominated by highly compact low-energy conformations which are entropically sup- 
pressed. Most promising for these studies appear sophisticated chain-growth methods based 
on Rosenbluth sampling 23] combined with improved pruning-enrichment strategies III, [28 ] 
which, in their original formulation, are particularly useful for the sampling in the G regime. 
For the analysis of the free zing transition, we apply in our simulations more generalized 



contact- density variants 



29 



301 ]. which have proven to be very successful in the low-energy 



sampling of protein 



solid substrates 



31 



ike heteropolymers 



29l | and the adsorption of polymers and peptides to 



32j. The precision of these algorithms when applied to lattice polymers 



as in the present study, is manifested by unraveling even finite-length effects induced by 
symmetries of the underlying lattice. 

The rest of the paper is organized as follows. In Sect. [Til the lattice model for flexible 
polymers and the employed chain-growth methods, which allow for a precise statistical 
sampling even in the low-temperature regime, are described. The conformational transitions 
the polymers experience on sc and fee lattices are discussed in Sect. IIHI Here, we also present 
our results for the scaling of the collapse transition temperature in comparison with various 
approaches known from the literature. Eventually, in Sect. HVl the paper is concluded by a 
summary of our findings. 



II. MODEL AND METHODS 



For our studies, we employ the interacting self-avoiding walk (ISAW) model for lattice 
polymers. In this model, the polymer chain is not allowed to cross itself, i.e., a lattice site 
can only be occupied by a single monomer. In order to mimic the "poor solvent" behavior 
in the energetic regime, i.e., at low temperatures, nearest-neighbor contacts of nonadjacent 
monomers reduce the energy. Thus, the most compact conformations possess the lowest 
energy. Formally, the total energy of a conformation X = {xi,X2, . . . , x^} of a chain with 
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N beads is simply given as 

E(X) = -e n NN (X), (1) 

where Eq is an unimportant energy scale (which is set 6q = 1 in the following) and wnn(X) 
is the number of nearest-neighbor contacts between nonbonded monomers. 

The total number of self-avoiding lattice conformations with m = N — 1 bonds scales 



as 



33, 



m 

C w ~/i w m 7 -\ (2) 

where /i is the effective coordination number of the lattice and 7 ~ 1.16 a universal ex- 
ponent. For the sc lattice, the connectivity constant is fi sc ~ 4.684 and in the fee case 



/if cc ~ 10.036 [35|, |36|, |37J, [38|, |39|, |40j, |4lj. Due to this exponential growth of the number of 



conformations, the investigation of all conformational transitions a homopolymer of a given 
length can experience requires employing numerical methods being capable of estimating 
the density of states for all possible energies with high precision. There exist mainly two 
strategies for generating and updating conformations within the stochastic search schemes 
typically used. Applying standard Markov chain Monte Carlo methods, conformational up- 
dates include semilocal changes of bond orientations as, among others, corner and end flips, 
crankshaft moves, and more non-local updates such as pivot rotations. In our work, we have 
used the alternative concept of chain growth. Depending on the lattice constraints, a new 
monomer is tried to be attached to an end of the already existing chain until the total length 
or a "dead end" is reached, i.e., all neighbors are already occupied and, thus, the chain end 



is trapped. In an early approach, the Rosenbluth method [271 ] . firstly the number of free 
nearest neighbors k for the possible placements of the Zth monomer is determined. Then, 
one of the possibilities (if any) is selected randomly. The peculiarity is that this algorithm 
introduces a bias as the (athermal) probability of generating a certain chain conformation 
X of length N, p(X) = Yl^^f 1 ' depends on the growth direction. Thus, identical walks 
can possess different construction probabilities, if they, e.g., were grown "forward" or "back- 
ward" . For correct statistics, this bias must be corrected by introducing Rosenbluth weights 
to(X) = p _1 (X). Actually, this bias can be utilized to increase the efficiency of the method 
in generating self-avoiding walks which is particularly useful for the ISAW model of lattice 
polymers. For this purpose it is convenient to introduce for each chain a thermal Rosenbluth 
weight W(X) = YI1L2 hi exp{ — [E{xj) — £ , (x i _ 1 )]//c B T} (where ks is the Boltzmann constant 
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which we set to ks = 1 in our analysis). Since this method is a clever kind of simple sam- 
pling, the partition sum can be estimated absolutely as Zn ~ J2i=i W(Xi)/M, where In is 
the number of successfully generated polymer chains of length N in M growth starts 11]. 
In principle, since Z^ = Y2e 9n{E) exp(— E/ksT), an absolute estimate of the density of 
states gjv(E), i.e., the degeneracy of energetic states E, is then also known. This is partic- 
ularly important for heteropolymers, where the ground-state degeneracy is considered as a 
measure for the stability of native folds of lattice proteins [29( . An essential improvement of 
the efficiency of this chain-growth approach was reached by the introduction of the pruned- 
enriched Rosenbluth method (PERM) Hj and its enhanced variants 28J], which combine 
Rosenbluth chain growth with a 'Go with the Winners' 42| strategy. In these algorithms, 
at each stage of the growth process copies of the already existing chain segment are created 
and continue growing independently, if the accumulated Rosenbluth weight is larger than 
an upper threshold. If the weight falls below a lower bound, the chain is pruned with a cer- 
tain probability (which is typically 1/2). Otherwise, the growth of the single chain si mply 



continues. This method has frequently been applied in studies of the point 
In our study, we use the nPERMss (new PERM with simple sampling) variant [28] for the 
simulation of G polymers with chain lengths of up to 32 000 (sc) and 4 000 (fee) monomers, 
respectively. 

For the analysis of the conformational behavior below the G point, we use even more 
sophisticated, generalized-ensemble variants which are independent of the temperature and 
yield an improved estimate for the density of states g(E) within a single simulation. These 
algorithms combine PERM-based chain growth with multicanonical 29( or flat-histogram 
techniques [30j and increase, in particular, the sampling of entropically suppressed ( "rare" ) 
conformations, which are, for example, essential for the study of the freezing transition. 
Due to the much higher demands in this regime, maximum chain lengths, for which precise 
results were reliably obtained, are N = 125 (sc) and 56 (fee), respectively. 

For our statistical analysis, it is convenient to define energetic statistical expectation 
values via the density of states, i.e., (0(E)) = J2e 9n(E)0(E) exp(—E /ksT) j 'Z^. The 
main results of our analysis are based on the peak structure of the specific heat which is 
defined as C V (T) = d(E)/dT = {(E 2 ) - (E) 2 ) /k B T 2 . 
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FIG. 1: Examples of specific-heat curves (per monomer) for a few exemplified short homopolymers 
on the sc lattice. Absolute errors (not shown) are smaller than 0.03 in the vicinity of the low- 
temperature peaks and smaller than 10 -5 in the onset of the O-transition region near T ~ 1.5. 

III. RESULTS AND DISCUSSION 

It was recently found for a bond-fluctuation model with inter-monomeric interaction ra- 
dius r = y/E that in the infinite chain-length limit collapse and freezing are indistinguishable 



phase transitions appearing at the same temperature (the G temperature Tq) 12|, [l3|. In 
a bead-spring FENE model analysis [14J, this phenomenon could not be observed: Both 
transitions exist in the thermodynamic limit and the crossover peaks in the specific heat 
remain well-separated. The same observation was made independently in Ref. [13] for the 
bond-fluctuation model with increased interaction range. In the following, we perform for 
the lattice polymer model (OQ) a detailed analysis of these transitions on regular sc and fee 
lattices and discuss the expected behavior in the thermodynamic limit. 



A. The expected peak structure of the specific heat 

Statistical fluctuations of the energy, as expressed by the specific heat, can signalize ther- 
modynamic activity. Peaks of the specific heat as a function of temperature are indicators 
for transitions or crossovers between physically different macrostates of the system. In the 
thermodynamic limit, the collective activity, which influences typically most of the system 
particles, corresponds to thermodynamic phase transitions. For a flexible polymer, three 
main phases are expected: The random-coil phase for temperatures T > Tq, where con- 
formations are unstructured and dissolved; the globular phase in the temperature interval 



6 



£=-79 E = -6Q 




E = 



FIG. 2: Representative conformations of a 64-mer in the different pseudophases: (a) Excitation 
from the perfect 4x4x4 cubic ground state (not shown, E = —81) to the first excited crystal 
state, (b) transition towards globular states, and (c) dissolution into random-coil conformations. 

T m < T < T@ (T m : melting temperature) with condensed, but unstructured ("liquid") 
conformations dominating; and for T < T m the "solid" phase characterized by locally crys- 
talline or amorphous metastable structures. In computer simulations, only polymers of finite 
length are accessible and, therefore, the specific heat possesses typically a less pronounced 
peak structure, as finite-length effects can induce additional signals of structural activity 
and shift the transition temperatures. These effects, which are typically connected with 
surface-reducing monomer rearrangements, are even amplified by steric constraints in lat- 
tice models as used in our study. Although these pseudotransitions are undesired in the 
analysis of the thermodynamic transitions, their importance in realistic systems is currently 
increasing with the high-resolution equipment available in experiment and technology. The 
miniaturization of electronic circuits on polymer basis and possible nanosensory applications 
in biomedicine will, therefore, require a more emphasized analysis of the finite-length effects 
in the future. 
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FIG. 3: Map of specific-heat maxima for several chain lengths taken from the interval N £ [8, 125]. 
Circles (0) symbolize the peaks (if any) identified as signals of the collapse (Tc-max > 1). The 
low-temperature peaks (+) belong to the excitation/freezing transitions (Tc-max < 0.8). The group 
of points in the lower left corner corresponds to polymers with N c + 1 monomers, where N c denotes 
the "magic" lengths allowing for cubic or cuboid ground-state conformations (see Fig. [Hand text). 

B. Simple-cubic lattice polymers 

Figure Q] shows typical examples of specific heats for very short chains on the sc lattice 
and documents the difficulty of identifying the phase structure of flexible homopolymers. 
The 27-mer exhibits only a single dominating peak - which is actually only an sc lattice 
effect. The reason is that the ground states are cubic (3x3x3) and the energy gap towards 
the first excited states is AE = 2 j^. Actually, also the most pronounced peaks for N = 48 
(4x4x3) and N = 64 (4 x 4 x 4) are due to the excitation of perfectly cuboid and cubic 
ground states, respectively. The first significant onset of the collapse transition is seen for 
the 48-mer close to T ~ 1.4. A clear discrimination between the excitation and the melting 
transition is virtually impossible. In these examples, solely for N = 64 three separate peaks 
are present. The plots in Figs. E^a)d^c) show representative conformations in the different 
pseudophases of the 64-mer. Due to the energy gap, the excitations of the cubic ground 
state with energy E = —81 (not shown) to conformations with E = —79 [Fig. 12(a)] result 
in a pseudotransition which is represented by the first specific-heat peak in Fig. [U The 
second less-pronounced peak in Fig. [T] around T pa 0.6 — 0.7 signalizes the melting into 
globular structures, whereas at still higher temperatures T pa 1.5 the well-known collapse 
peak indicates the dissolution into the random-coil phase. 
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FIG. 4: (a) Collapse (0) and crystallization/excitation (+) peak temperatures of the specific heat 
for all chain lengths in the interval N € [8, 125], (b) values of the specific-heat maxima in the same 
interval. Error bars for the collapse transition data (not shown) are much smaller than the symbol 
size. peaks appear starting from iV = 41. For the sake of clarity, not all intermediate O data 
points are shown (only for N = 41, 45, 50, . . .). 

The distribution of the maximum values of the specific heat C™ ax with respect to the 
maximum temperatures Tc-max is shown in Fig. Not surprisingly, the peaks belonging 
to the excitation and freezing transitions (+) appear to be irregularly "scattered" in the 
low-temperature interval < T^max < 0.8. The height of the peaks indicating the collapse 
transition of the finite-length polymers (©) is, on the other hand, monotonously increasing 
with the collapse-peak temperature. 

Figures H](a) andlU^b), showing the respective chain-length dependence of the maximum 
temperatures and maximum specific-heat values, reveal a more systematic picture. At least 
from the results for the short chains shown, general scaling properties for the freezing transi- 
tion cannot be read off at all. The reason is that the low-temperature behavior of these short 
chains is widely governed by lattice effects. This is clearly seen by the "sawtooth" segments. 
Whenever the sc chain possesses a "magic" length N c such that the ground state is cubic or 
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FIG. 5: Examples of specific heats for polymers with (a) N c = 4 x 3 x 3 = 36 and N c ± 1 and 
(c) JV c = 4x4x4 = 64 and N c ± 1 monomers. In (b) and (d), the respective densities of states 
are shown (lines are only guides to the eye). Symbols *, + and x emphasize the lowest-energy 
states. Note the energy gaps between the ground and the first excited states in the compact cases 
N c = 36, 64 and the dip for N = 35, 63. 

cuboid (i.e., N c e Af c = {8, 12, 18, 27, 36, 48, 64, 80, 100, 125, ...}), the energy gap AE = 2 
between the ground-state conformation and the first excited state entails a virtual energetic 
barrier which results in an excitation transition. Since entropy is less relevant in this regime, 
this energetic effect is not "averaged out" and, therefore, causes a pronounced peak in the 
specific heat [see Fig. IU(b)] at comparatively low temperatures [Fig. H(a)]. This peculiar sc 
lattice effect vanishes widely by increasing the length by unity, i.e., for chain lengths iV c + 1. 
In this case, the excitation peak either vanishes or remains as a remnant of less thermody- 
namic significance. The latter appears particularly in those cases, where N = N c + 1 with 
N c = L 3 (with L being any positive integer) is a chain length allowing for perfectly cubic 
ground states. Increasing the polymer length further, the freezing peak dominates at low 
temperatures. Its peak increases with the chain length, whereas the peak temperature de- 
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creases. Actually, with increasing chain length, the character of the transition converts from 
freezing to excitation, i.e., the entropic freedom that still accompanies the melting/freezing 
process decreases with increasing chain length. In other words, cooperativity is lost: only a 
small fraction of monomers - residing in the surface hull - is entropically sufficiently flexible 
to compete the energetic gain of highly compact conformations. This flexibility is reduced 
the more, the closer the chain length TV approaches a number in the "magic" set M c - If the 
next length belonging to M c is reached, the next discontinuity in the monotonic behavior 
occurs. Since noticeable "jumps" are only present for chain lengths whose ground states 
are close to cubes (N c = L 3 ) or cuboids with N c = L 2 (L ± 1), the length of the branches 
in between scales with AiV c ~ L 2 ~ N^ 3 . Therefore, only for very long chains on the sc 
lattice, for which, however, a precise analysis of the low-temperature behavior is extremely 
difficult, a reasonable scaling analysis for Tc^(N) and C™ ax (iV) could be performed. 

Exemplified for chains of lengths N = N c — 1, N c , and N c + 1 with N c = 36, 64, specific 
heats and densities of states are shown in Figs.[5]^a)-[5](d), which exhibit the length-dependent 
characteristic properties discussed above. While in Fig. [5](a) for chain lengths around N = 36 
only some low-temperature activity is visible and the collapse transition is vaguely indicated 
by a broad shoulder, the transition characteristics are better resolved for N = 64 shown in 
Fig. EJ^c). The most pronounced low-temperature peak of the specific heat of the 64-mer, for 
example, is the excitation peak, the second peak belongs to the freezing transition, and the 
third, still very shallow peak signals the collapse transition. The low-temperature behaviors 
of the 63-mer and the 65-mer are quite different: While the low-temperature peak of the 63- 
mer close to T pa 0.4 is due to excitation as is indicated by the E = —79 "dip" in the density 
of states in Fig. [5]^d) [similar to the 35-mer in Fig. [5^b) at E = —38], the relevant peak for 
the 65-mer is the freezing peak close to T pa 0.6. In this case, the excitation is of much less 
relevance (although it is still reflected by a small peak near T pa 0.3). This is a consequence 
of the missing convex lowest-energy dip in the density of states (or microcanonical entropy). 
The convex monotony is a signal of a strong first-order phase separation 44j, 145]. This is 
confirmed by analyzing the canonical energy distributions for the examples N = N c — 1, N c , 
and N c + 1 with N c = 36, 64 shown in Figs. [6](a)-[6](c) for temperatures close to the respective 
excitation and freezing transitions. For N c — 1 = 35, 63 [Fig. EJ^a)], the pronounced excitation 
transition is expressed by the respective double peaks with the strong gap in between, which 
are for the polymers with chain lengths N c = 36, 64 [Fig. MJo)} due to the energy gap 



11 



between ground state and first excited state. This induces the first-order-like character of 
this pseudotransition. The energy distributions for various temperatures shown in Fig. [6](c) 
for the case N = N c + 1 = 65 do not exhibit, on the other hand, pronounced double-peaked 
shapes. The excitation transition at extremely low energies is still weakly present as a small 
shoulder in the distribution at the corresponding temperature. The freezing transition is 
associated with slightly larger energies (and temperature) and visible in the distribution 
with a weak tendency to a double-peaked shape. 

The collapse transitions of the finite-length polymers are not affected by the intricate 
low-energy conformations on the sc lattice and exhibit a continuous monotony. This will be 
analyzed in Sect. IIIIDl in more detail. 

C. Polymers on the fee lattice 

The general behavior of polymers on the fee lattice is comparable to what we found 
for the sc polymers. The main difference is that excitations play only a minor role, and 
the freezing transition dominates the conformational behavior of the fee polymers at low 
temperatures. Nonetheless, finite-length effects are still apparent as can be seen in the 
chain-length dependence of the peak temperatures and peak values of the specific heats 
plotted in Fig. [7](a) and Fig. [7(b), respectively. Figure [7](a) shows that the locations of the 
freezing and collapse transitions clearly deviate with increasing chain lengths and we hence 
can conclude that also for fee polymers there is no obvious indication that freezing and 
collapse could fall together in the thermodynamic limit. 

Similar to the sc polymers, the finite-length effects at very low temperatures are ap- 
parently caused by the usual compromise between maximum compactness, i.e., maximum 
number of energetic (nearest-neighbor) contacts, and steric constraints of the underlying 
rigid lattice. The effects are smaller than in the case of the sc lattice, as there are no obvi- 
ous "magic" topologies in the fee case. Ground-state conformations for a few small polymers 
on the fee lattice are shown in Fig. [HJ The general tendency is that the lowest-energy con- 
formations consist of layers of net planes with (111) orientation, i.e., the layers themselves 
possess triangular pattern with side lengths equal to the fee nearest-neighbor distance V2 
(in units of the lattice constant). This is not surprising, as these conformations are tightly 
packed which ensures a maximum number of nearest-neighbor contacts and, therefore, lowest 
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FIG. 6: Energy distributions P(E) at low temperatures for sc lattice polymers with (a) N = 
N c — 1 =35 and 63, (b) N = N c = 36 and 64 monomers. In (c), P(E) is shown for the chain 
with length N = N c + 1 = 65 for several temperatures T = 0.2, 0.3, . . . , 0.9. The distributions for 
iV = 37 (not shown) are similar. Note that lines are only guides to the eyes. 

conformational energy. An obvious example is the ground-state conformation of the 13-mer 
as shown in Fig. 0(a) which corresponds to the intuitive guess for the most closely packed 
structure on an fee lattice: a monomer with its 12 nearest neighbors ("3-7-3" layer struc- 
ture). A simple contact counting yields 36 nearest-neighbor contacts which, by subtracting 
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FIG. 7: Peak temperatures (a) and peak values (b) of the specific heat for all chain lengths 
N = 8, . . . , 56 of polymers on the fee lattice. Circles (0) symbolize the collapse peaks and low- 
temperature peaks (+) signalize the excitation/freezing transitions. The error bars for the collapse 
transition are typically much smaller than the symbol size. Only for the freezing transition of 
longer chains, the statistical uncertainties are a little bit larger and visible in the plots. peaks 
appear starting from N = 19. For clarity, @ data points are only shown for N = 19, 25, 30, .... 

the N — l = 12 covalent (nonenergetic) bonds, is equivalent to an energy E = —24. However, 
this lowest-energy conformation is degenerate. There is another conformation (not shown) 
consisting of only two "layers", one containing 6 (a triangle) and the other 7 (a hexagon) 
monomers ("6-7" structure), with the same number of contacts. 

A special case is the 18-mer. As Fig. WJd) shows, its ground state is formed by a complete 
triangle with 6 monomers, a hexagon in the intermediate layer with 7 monomers, and an 
incomplete triangle (possessing a "hole" at a corner) with 5 monomers ("6-7-5" structure). 
Although this imperfection seems to destroy all rotational symmetries, it is compensated 
by an additional symmetry: Exchanging any of the triangle corners with the hole does not 
change the conformation at all! Thus, the seeming imperfection has a similar effect as the 
energetic excitation and causes a trivial entropic transition. This explains, at least partly, 
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iV = 13 iV = 18 iV = 19 

Eo 24 ^ 39 Eq = -42 

iV = 27 
£ij = -67 

FIG. 8: Ground-state conformations and energies of the (a) 13-, (b) 18-, (c) 19-, (d) 27-, and (e) 
30-mer on the fcc-lattice (bonds not shown). 

why the 18-mer exclusively exhibits an additional peak in the specific heat at very low 
temperatures [see Fig. Et^a)]. A similar reasoning presumably also applies to the anoma- 
lous low-temperature peaks of the 32-, 46-, and 56-mers, but for these larger ground-state 
conformations it does not make much sense to go into such intricate details. 

The expectation that the 19-mer, which can form a perfect shape without any "holes" 
("6-7-6" structure), is a prototype of peculiar behavior is, however, wrong. This is due to 
the existence of degenerate less symmetric ground-state conformations [as the exemplified 
conformation in Fig. [H^c)]. 

The described geometric peculiarities are, however, only properties of very short chains. 
One of the largest of the "small" chains that still possesses a non-spherical ground state, is 
the 27-mer with the ground-state conformation shown in Fig. E^d). For larger systems, the 
relative importance of the interior monomers will increase, because of the larger number of 
possible contacts. This requires the number of surface monomers to be as small as possible 
which results in compact, sphere- like shapes. A representative example is the 30-mer shown 
in Fig.EKe). 
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TABLE I: Tq values on the sc and fee lattice from literature. 



lattice type 








model 


Ref. 


sc 


3.64 




4.13 


single chain 


[201 




3.713 


± 


0.007 


single chain 


[211 




3.650 




0.08 


single chain 


[231 




3.716 


± 


0.007 


single chain 


[8] 




3.60 


± 


0.05" 


single chain 


[9] 




3.62 


± 


0.08 6 


single chain 


[101 




3.717 


± 


0.003 


single chain 


[111 




3.717 


± 


0.002 


polymer solution 


[161 




3.745 






lattice theory 


[25] 




3.71 


± 


0.01 


polymer solution 




fee 


8.06 




9.43 


single chain 


[20] 




8.20 


± 


0.02 


single chain 


[22, 24J 




8.264 






lattice theory 


[25] 



"Origina 
b In Ref. 



ilyigiven as (3 e = 0.2779 ± 0.0041 |9| 
10J given as (3 e = 0.276 ± 0.006. 



D. The transition revisited 



The scaling behavior of several quantities at and close to the G point in three dimensions 
has been the subject of a large number of field-theoretic and computational studies 0, 



11 



16 



17 



18 



21 



22 



23 



24 



25 



261 ] . Nonetheless, the somewhat annoying result is 
that the nature of this phase transition is not yet completely understood. The associated 
tricritical lim„^ 0(n) field theory has an upper critical dimension d c = 3, but the predicted 
logarithmic corrections [4, f], [f| could not yet be clearly confirmed from the numerical data 
produced so far. In our study of freezing and collapse on regular lattices, we mainly focused 
on the critical temperature Tq for polymers on the sc and on the fee lattice. The sc value of 
Te has already been precisely estimated in several studies, but only a few values are known 
for the fee case. Some previous estimates in the literature are compiled in Table [B 

As our main interest is devoted to the expected difference of the collapse and freezing 
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temperatures, we will focus here on the scaling behavior of the finite-size deviation of the 
maximum specific-heat temperature of a finite- length polymer from the 6 temperature, 
T C (N) — Tq, as it has also been studied for the bond-fluctuation model 



off-lattice FENE polymer 14|, as well as for polymer solution models 
latter case, Flory-Huggins mean-field theory [46( suggests 
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2. 



13, 



l 



l 



l 



l 

2iV' 



-3| and the 
In the 



(3) 



T clit (N) T v/jV 

where T cr i t (iV) is the critical temperature of a solution of chains of finite length N and 
Tq = liniTv^oo Trait (iV) is the collapse transition temperature. In this case, field theory ^] 
predicts a multiplicative logarithmic correction of the form T cr ; t (A^) — Te ~ iV -1 / 2 [ln iV] -3 / 11 . 
Logarithmic corrections to the mean-field theory of single chains are known, for example, 
for the finite-chain Boyle temperature Tg(iV), where the second virial coefficient vanishes. 
The scaling of the deviation of Tb(N) from Tq reads js]: 

1 



T B (N)-T ( 



e 



(4) 



^(ln iV) 7 /n 

n 

In Ref. [14|, it is claimed that, for their data obtained from simulations with the FENE 
potential, this expression can also be used as a fit ansatz for T C (N) — Tq. However, also the 
mean-field-motivated fit without explicit logarithmic corrections, 

a 2 



T C (N)-T 



e 



ai 



(5) 



has been found to be consistent with the off-lattice data [ijj], and also with the re- 
sults obtained by means of the bond- fluctuation model of single chains with up to 512 



monomers 



12 



131 ] . Up to corrections of order N 3//2 , Eq. (jSJ) is equivalent to 
1 1 Si 



T C (N) 



1 



02 

N + N' 



(6) 



which was found to be consistent with numerical data obtained in grandcanonical analyses 
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13, 



it 



of lattice homopolymers and the bond-fluctuation model 

The situation remains diffuse as there is still no striking evidence for the predicted loga- 
rithmic corrections (i.e., for the field-theoretical tricritical interpretation of the O point) from 
experimental or numerical data. Using our data from independent long-chain nPERMss 28 ] 
chain-growth simulations (sc: iV max = 32 000, fee: iV max = 4 000) in the vicinity of the col- 
lapse transition, we have performed a scaling analysis of the iV-dependent collapse transition 
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0.015 0.03 0.045 0.06 0.075 0.09 0.105 



1/ViV 



FIG. 9: Inverse collapse temperatures for several chain lengths on sc (N < 32 000) and fee lattices 
(N < 4 000) . Drawn lines are fits according to Eq. ([6]) . 

temperatures T C (N), identified as the collapse peak temperatures of the individual specific- 
heat curves, and estimated from it the N — > oo limit Tq. For the single-chain system, field 
theory predicts the specific heat to scale at the 9 point like CV(T = Tq)/N ~ (In iV) 3 / 11 . 



Short-chain simulations 23] did not reveal a logarithmic behavior at all, whereas for long 
chains a scaling closer to In N was read off [llj. The situation is similar for structural 
quantities such as the end-to-end distance and the gyration radius. Figure [9] shows our data 
points of the inverse collapse temperature T c _1 from the simulations on the sc (left scale) 
and on the fee lattice (right scale), plotted against iV" 1 / 2 . Error bars for the individual data 
points in Fig. [9] were obtained by jackknife error estimation |47J] from several independent 
simulation runs. Also shown are respective fits according to the ansatz ([6]). Optimal fit 
parameters using the data in the intervals 200 < N < 32 000 (sc) and 100 < iV < 4 000 (fee) 
were found to be = 3.72(1), 6i « 2.5, and a 2 « 8.0 (sc) and T| c = 8.18(2), a x « 1.0, 
and a2 ~ 5.5 (fee). In addition, we investigated also other fit functions motivated by field 
theory and mean- field-like approaches, corresponding to Eqs. (EI) (IHJ) , each of which also 
with different fit ranges. These results are listed in Tables [Til and respectively. In order 
to decide which of the fits is consistent with our data, the x 2 test is used. Depending on 
the sizes of the data sets entering into the analyses and the number of fit parameters, there 
are 2 to 6 degrees of freedom df. We make the typical assumption that deviations of the fit 
from the used data set are significant, if x 2 > xl f; o.o5' i- e -' ^ X 2 ues i n the 5% tail of the 
Pd { (x 2 ) distribution of x 2 values. In this case, with 95% probability the deviations between 
data and fit function are not random. The thresholds for the different degrees of freedom 
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TABLE II: Values of Tq on the sc lattice from different fits and their x 2 tests with several degrees 
of freedom d{. 



fit function 




To 


a resp. oti, 0,2 




A 


Y 2 /dr 


df 

1 






3 6671 ± 0051 


—26.0 


500 


- 32 000 


12.7 


4 


TAN) T© = ■ 


a 


3 6741 ± 0053 


—26.5 


1 000 


- 32 000 


6.61 


3 






3.6898 ± 0.0065 


-28.3 


2 000 


- 32 000 


1.39 


2 






3 7353 ± 0056 


—8.24 


500 


- 32 000 


0.57 


4 




a fir, N\ 7 / n 


3.7370 ± 0.0057 


-8.27 


1 000 


- 32 000 


0.21 


3 






3.7398 ± 0.0072 


-8.36 


2 000 


- 32 000 


0.10 


2 






3 6164 ± 0048 


—83.0 


500 


- 32 000 


39.5 


4 




a 


3.6287 ± 0.0049 


-86.0 


1 000 


- 32 000 


20.7 




/A(lnA) 7 /n 


3 






3.6531 ± 0.0059 


-97.6 


2 000 


- 32 000 


4.14 


2 






3 7255 ± 0060 

O • 1 ^— '.) 'J _1_ \Ji\J\J\J\J 


2.6 


500 


- 32 000 


0.28 


4 


1 1 


4 1 + M 


q 7945 + n f)061 


2 6 


1 000 


- 32 000 


1 1 


'J 


T C (N) T 


\y/N 2Nj 






3.7221 ± 0.0076 


2.6 


2 000 


- 32 000 


0.03 


2 






3.7173 ±0.0071 


2.5, 8.0 


200 


- 32 000 


0.05 


4 


1 1 




3.7173 ±0.0104 


2.5, 8.0 


500 


- 32 000 


0.07 


3 


T C (N) T 


v^V 2A 






3.7194 ±0.0131 


2.5, 6.3 


1000 


- 32 000 


0.07 


2 






3.7030 ± 0.0059 


-32, 135 


200 


- 32 000 


0.53 


4 


T C (N)-T e = - 


ai a 2 
/A A 


3.7090 ± 0.0078 


-32, 161 


500 


- 32 000 


0.25 


3 






3.7140 ±0.0104 


-33, 186 


1000 


- 32 000 


0.12 


2 



lie between xi-o.os/^f = 3-0 and X&o.os/di = 2-1- The calculated x 2 values associated with 
the data sets and the fit functions used are also listed in Tables [TT] and IIHI 

From the results in Table [III for the polymers on the sc lattice, we find that the two- 
parameter mean-field- like fits (|SJ) and ([6]) as well as the single-parameter fit according to (J3j) 
are consistent with our data. Surprisingly poor, on the other hand, is the goodness of 
the fit against the logarithmic scaling (J3J). Even more astonishing is, however, the good 
coincidence with a logarithmic fit of the "wrong " form A" -1 / 2 (In AO 7/11 with the data. Sum- 
marizing these results, if logarithmic corrections as predicted by tricritical field theory are 
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TABLE III: Values of T@ on the fee lattice using the same methodology as in Table [II] for the sc 
lattice. 



fit function 




T e 


a resp. a\, a 2 




N 


x 2 M 


<k 






7.2673 ± 0.0052 


-34.1 


100 


- 4000 


1000 


6 






7.5760 ± 0.0070 


-39.4 


150 


- 4 000 


418 


5 


T C (N) - T© = ■ 


a 
\/N 


7.7101 ± 0.0080 


-42.0 


210 


- 4000 


213 


4 






7.8445 ± 0.0096 


-45.5 


300 


- 4 000 


67.8 


3 






7.9561 ± 0.0013 


-50.0 


500 


- 4000 


13.6 


2 






7.9218 ± 0.0064 


-15.2 


100 


- 4000 


200 


6 






8.0757 ± 0.0083 


-16.1 


150 


- 4000 


69.2 


5 


T C (N)-T e = - 


n 

" (In AO 7/11 


8.1356 ±0.0093 


-16.5 


210 


- 4000 


32.7 


4 






8.1953 ±0.0110 


-17.0 


300 


- 4 000 


9.45 


3 






8.2468 ± 0.0149 


-17.6 


500 


- 4 000 


1.00 


2 






6.8260 ± 0.0045 


-79.5 


100 


- 4 000 


2 000 


6 








— QQ 3 


lull 


^± uuu 


i nnn 


K 
o 


T C {N)-T e = - 


a 


7.4166 ± 0.0072 


-110.2 


210 


- 4 000 


500 


4 


s/NilnN) 7 / 11 






7.6051 ± 0.0087 


-125.7 


300 


- 4000 


164 


3 






7.7544 ± 0.0011 


-146.3 


500 


- 4 000 


38.1 


2 






8.5434 ± 0.0110 


1.27 


100 


- 4 000 


111 


6 






8.4208 ± 0.0120 


1.23 


150 


- 4000 


29.6 


5 


1 1 


a \ /AT + 9/V ) 


8.3821 ±0.0125 


1.22 


210 


- 4 000 


13.5 


4 


ic{w ) 


V V 2iV / 






8.3369 ±0.0141 


1.20 


300 


- 4 000 


3.00 


3 






8.3048 ±0.0187 


1.18 


500 


- 4000 


1.15 


2 






8.1778 ±0.0169 


1.04, 5.49 


100 


- 4 000 


0.81 


5 


1 1 


ai a2 


8.1987 ±0.0211 


1.06, 5.04 


150 


- 4 000 


0.32 


4 


T C (N) T 


^A 7 + 2N 


o.zwi ± u.uzoy 


l.U / , 4. ( 




a nnn 

— 4 UUU 


n 91 


q 






8.2288 ± 0.0386 


1.09, 4.18 


300 


- 4 000 


0.11 


2 






8.0374 ± 0.0110 


-58.9, 360 


100 


- 4 000 


13.1 


5 


T C (N)-T & = - 


a\ a 2 


8.0876 ± 0.0133 
8.1219 ±0.0163 


-61.5, 414 
-63.5, 461 


150 
210 


- 4 000 

- 4 000 


4.71 
1.81 


4 
3 






8.1640 ±0.0244 


-66.5, 541 


300 


- 4000 


0.04 


2 
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present at all, even chain lengths N = 32 000 on an sc lattice are too small to observe de- 
viations from the mean-field picture. At least, the goodness of the logarithmic fit with the 
"wrong" exponent +7/11 could lead to the speculative conclusion that for iV < 32 000 mul- 
tiplicative and additive logarithmic corrections to scaling are hidden in the fit parameters 
of the "mean-field-like fits" . The subleading additive corrections are expected to be of the 
form In (In N)/\n 2 N [3]. They thus not only disappear very slowly - they are also even of 
the same size as the leading scaling behavior, which makes it extremely unlikely to observe 
the logarithmic corrections in computational studies at all Similar additive logarithmic 
scaling is also known, for example, from studies of the two-dimensional XY spin model 48]. 
The estimated sc temperatures from the good fits are in perfect agreement with the most 
reliable estimates from literature. 

The corresponding fee results are listed in Table IIHI In this case, only the fit function ([6]) 
is independent of the data sets used and, therefore, consistent with the data obtained for 
all chain lengths. However, the noticeable improvement of the goodness for the fits to (jHJ), 
([3]) and the "wrong" A~ 1//2 (ln A) 7 / 11 form by excluding the very short chains from the data 
sets considered, leads to the conclusion that even chains with A = 4 000 monomers on the 
fee lattice are also too short to find evidence for the logarithmic corrections to mean-field 
scaling. Our best estimates for the fee temperature agree nicely with the results from 



Refs. 



22 



24] . 



IV. SUMMARY 



Employing sophisticated chain-growth algorithms, we have performed computer simula- 
tions of homopolymers on sc and fee lattices in order to analyze freezing and collapse of these 
chains. Particular attention has been devoted to the question whether these transitions fall 
together in the thermodynamic limit as it was reported recently from similar studies of a 
specific bond-fluctuation model. In our analysis, we focus on the shifts of the specific-heat 
peaks in dependence of the chain-lengths considered. 

For polymers on the sc lattice, we find a remarkably systematic pattern of the freezing 
transition which can be explained by lattice effects of the finite-length systems. In fact, 
the high precision of our data allows us to reveal a noticeable difference in the behavior 
of "magic" chain lengths that allow for cubic or cuboid conformations. In these cases, an 
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energy gap exists between the ground-state conformations and the first excitations. This 
peculiarity causes a first-order- like pseudotransition which is typically more pronounced 
than the separate freezing transition. Surprisingly, this effect vanishes widely for polymers 
with slightly longer chain lengths. The freezing temperature decreases with increasing chain 
length until the next "magic" length is reached. Polymers on the fee lattice behave similarly, 
but the relevant geometries are more complex. 

We have also performed comprising analyses of the collapse temperature deviations for 
finite-length polymers from the G temperature, i.e., the collapse transition temperature 
in the limit of infinite chain-length. We studied chains with lengths of up to 32 000 (sc) 
and 4 000 (fee) monomers, respectively. The thus obtained data were fitted against several 
fit functions motivated by field-theoretic and mean-field-like approaches. For the chain 
lengths studied, we find no evidence for logarithmic corrections as predicted by tricritical 
field theory. We conclude that the chain lengths are still too short to uniquely identify 
logarithmic corrections which are probably effectively taken into account by the amplitudes 
of the dominant mean-field terms. 

From our results for the freezing and the collapse transition, we conclude that both 
transitions remain well separated also in the extrapolation towards the thermodynamic limit. 
This is the expected behavior as it is a consequence of the extremely short range of attraction 
in the nearest-neighbor lattice models used. Considering a more general square-well contact 
potential between nonbonded monomers in our parametrization, 



oo r < 1, 



v(r) 



-1 1 < r < A, (7) 
A < r, 



the attractive interaction range is simply R — A — 1. In our single-chain study of sc and 
fee lattice models, we have A — ^ 1 and thus R — > 0. Since this R value is well below 
a crossover threshold known for colloids interacting via Lennard- Jones-like and Yukawa 
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50, 



5l|, we interpret our 



potentials, where different solid phases can coexist, R { c ] « 0.01 
low-temperature transition as the restructuring or "freezing" of compact globular shapes 
into the (widely amorphous) polymer crystals. 

Following Ref. 52|, there is also another phase boundary, namely between stable and 
metastable colloidal vapor-liquid (or coil-globule) transitions, in the range 0.13 < Rr < 
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0.15. Other theoretical and experimental approaches yield slightly larger values, 



0.25 



49 



53 



54 



55]. Below Fi¥\ the liquid (globule) phase is only metastable. The specific 
bond-fluctuation model used in Ref. [3] corresponds to R = 0.225, i.e., it lies in the crossover 
regime between the stable and metastable liquid phase [13). Consequently, the crystallization 
and collapse transition merge in the infinite-chain limit and a stable liquid phase was only 
found in a subsequent study of a bond-fluctuation model with larger interaction range [131 ]. 

Qualitatively, analogous to the behavior of colloids, our considerations would explain 
the separate stable crystal, globule, and random-coil (pseudo)phases that we have clearly 
identified in our lattice polymer study. Since the range of interactions seems to play a 
crucial, quantitative role, it is an interesting, still widely open question to what extent the 
colloidal picture in the compact crystalline and globular phases is systematically modified 
for polymers with different nonbonded interaction ranges, where steric constraints (through 
covalent bonds) are a priori not negligible. 
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